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■ Abstract 

Si 

' We develop a formalism for calculating cosmic microwave background (CMB) 

■ temperature and polarization anisotropies in cosmological models with Brans- 
T-H ■ Dicke gravity. We then modify publicly available Boltzmann codes to calcu- 

, late numerically the temperature and polarization power spectra. Results are 

\^ \ illustrated with a few representative models. Comparing with the general- 

' relativistic model with the same cosmological parameters, both the amplitude 

■ and the width of the acoustic peaks are different in the Brans-Dicke models. 

■ We use a covariance-matrix calculation to investigate whether the effects of 
Brans-Dicke gravity are degenerate with those of variation in other cosmolog- 
ical parameters and to simultaneously determine whether forthcoming CMB 

^ ^ maps might be able to distinguish Brans-Dicke and general-relativistic cosmol- 

' ogy. Although the predicted power spectra for plausible Brans-Dicke models 

c/3 . differ from those in general relativity only slightly, we find that MAP and/or 

• • , the Planck Surveyor may in principle provide a test of Brans-Dicke theory 

\ that is competitive to solar-system tests. For example, if all other parame- 

^ ' ters except for the CMB normalization are fixed, a value of the Brans-Dicke 

' parameter u; as large as 500 could be identified with MAP, and for Planck, 

values as large as u; ~ 3000 could be identified; these sensitivities are de- 
creased roughly by a factor of 3 if we marginalize over the baryon density, 
Hubble constant, spectral index, and reionization optical depth. In more gen- 
eral scalar-tensor theories, uj may evolve with time, and in this case, the CMB 
probe would be complementary to that from solar-system tests. 
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I. INTRODUCTION 



The Jordan-Fierz-Brans-Dicke theory (heretofore, we will call it Brans-Dicke the- 
ory for simplicity) is the simplest example of a scalar-tensor theory of gravity Q]. Recently, 
scalar-tensor theories have received renewed interest, because such theories are generic pre- 
dictions of superstring theory [Q and other higher-dimensional gravity theories 0. Fur- 
thermore, scalar-tensor theories have also found application in construction of inflationary 
models, including some models based on first-order phase transitions that evade the "graceful 
exit" problem [0-0. 

In Brans-Dicke theory, Newton's constant becomes a function of space and time, and 
a new parameter u is introduced. General relativity is recovered in the limit uj oo. 
Solar-system experiments using Viking ranging data [O have constrained u > 500. Recent 



VLBI measurement of the time delay of millisecond pulsars may further raise this limit [|12 
However, these experiments are all "weak-field" experiments and probe only a limited range 
of space and time. To effectively constrain more general scalar-tensor theories, one would 
also like to have "strong-field" experiments, such as that provided by the binary pulsar 



ri| , p^ . It was also pointed out [0] that in cosmological models based on more general 
scalar-tensor theories (in which u can vary), there is generally an attractor mechanism that 
drives a; to cx) at late times. Thus, it is possible that gravity differed considerably from 
general relativity in the early Universe, even if general relativity seems to work well today. 
Big-bang nucleosynthesis |TB|-P^ provides one test of gravity at early times. 

With the advent of precise new cosmic microwave background (CMB) data, it is natural 
to inquire whether the CMB might be able to provide new tests of Brans-Dicke theory (or of 
more general scalar-tensor theories). The Microwave Anisotropy Probe (MAP) [|1^] (to be 
launched in Fall 2000) and Planck Surveyor (to be launched in 2007) satellites as well as 



many ground-based and balloon-borne experiments will measure the CMB anisotropy with 
unprecedented precision, thus providing a wealth of information about the early Universe. 
The advantage of the CMB anisotropy is that it involves fairly simple linear physics and is 
thus very "clean" . Furthermore, the CMB anisotropy probes a different era of the cosmos. 
Thus, at least in principle one may see the presence of a scalar-tensor theory that has been 
driven by the attractor mechanism to the general-relativity limit in the current epoch |13 



The possibility of testing scalar-tensor gravity with CMB anisotropy has already been 
noted. For example, the original version of extended inflation was ruled out because the 
bubbles formed during the phase transition would have produced CMB anisotropy larger 
than that observed unless the Brans-Dicke parameter uj was less than 30 0. 

The general behavior of cosmological perturbations in Brans-Dicke cosmology was stud- 
ied analytically in Refs. PD|,pT| , but they did not consider realistic models. Peebles and Yu, 



in their pioneering study of CMB anisotropy |2^, considered a more realistic model with 



Brans-Dicke gravity, and they showed how the difference in the expansion rate affects the 



photon transfer function. More recently, Liddle et al. [|23| estimated that in Brans-Dicke 
theory, the epoch of radiation-matter equality is shifted, 
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and this accordingly affects the scale at which the present-day matter power spectrum turns 
over. 

In the particular case of cosmologies based on chaotic-inflation models, the production 
of fluctuations during inflation with scalar-tensor gravity has been studied in Refs. p^-p6|. 
They concluded that isocurvature perturbations could be produced during inflation, but are 
in general negligible compared with the adiabatic perturbations. In some inflation models, 
the spectrum of density perturbations may be affected, and for scalar-tensor theories with 
variable cu, the spectral index for primordial perturbations may change with scale. For 
example, in some Brans-Dicke inflationary models there is a slight tilt in the spectrum of 
density perturbations, and a limit on the variation of uj can be obtained from the COBE 



observation measurement of the spectral index [^, but only within the context of this very 
particular inflation model. 

In this paper, we perform a complete calculation of the CMB anisotropy in the Brans- 
Dicke theory. To do this we modify a standard code for CMB anisotropy calculation 
to accommodate Brans-Dicke theory. Our modified code may be used to calculate the 
anisotropy in any given cosmological model. Although our code can accommodate isocurva- 
ture perturbations as well, we present numerical results only for models with nearly scale- 
invariant spectra of primordial adiabatic perturbations for the following reasons: If acoustic 
peaks like those expected from adiabatic perturbations are observed, then it is plausible 
that we might understand structure formation well enough to use CMB anisotropy to look 
for tiny deviations from general relativity. If it appears that some more complicated physics 
gave rise to structure formation, then it is unlikely that the CMB will provide a precision 
tool to study gravity. 

We limit ourselves here to the simplest scalar-tensor theory: i.e., the original Brans-Dicke 
theory, for which the Brans-Dicke parameter uj is fixed. We will leave the more general case 
with variable uj to future work. Likewise, we concentrate on flat CDM models, including 
those with a cosmological constant, but without hot dark matter. 

This paper is organized as follows: In Section II, we develop the formalism for the 
calculation. In Section III, we study the behavior of the background cosmology and discuss 
the initial conditions for the perturbations in the Brans-Dicke field. Numerical results are 
presented in Section IV, and we also discuss the detectability of Brans-Dicke theory there. 
Section V then summarizes and concludes. We briefly describe the numerical implementation 
of the calculation in the Appendix. Throughout this paper, we use natural units, c = h = 



II. FORMALISM 
A. Brans-Dicke Theory 

The Lagrangian density for the Brans-Dicke theory is 



(2) 
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where $ is the Brans-Dicke field, and Lm is the Lagrangian density for the matter fields, 
whose equations of motion are not affected. For convenience, we also define a dimensionless 
field 

= (3) 

where G is the Newtonian gravitational constant measured today. 
The Einstein equations are generalized to 

G^. = ^T^. + - lgf.u^.,t) + ^{^■,^^^^ " 9^^u^^), (4) 

where T^jy is the stress tensor for all matter except for the Brans-Dicke field. The equation 
of motion for $ is 

m = — ^T. (5) 
Here T = is the trace of the energy-momentum tensor. 

B. Background Cosmology 

The unperturbed part of the metric in a fiat universe can be written as 

ds^ = a\T){-dT^ + 'jijdx'dx^). (6) 

where a is a function of the conformal time r only, and is the fiat-space metric. The 
unperturbed stress-energy tensor has components 

n = -p, T0 = 0, T)=]rf), (7) 

in the comoving frame. The equations describing the background evolution are 

p' + 3-(p + p) = 0, (8) 
a 

0" + 2^0' = |^(p-3p), (10) 

where the prime denotes derivative with respect to r, and p and p are the total density and 
pressure of the Universe, respectively. General relativity is recovered in the limits 

^0, $' ^ 0, u^oo. (11) 
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C. Cosmological Perturbations 



We can write the perturbed metric as 

9f,u = a^ilpiv + hf,^), (12) 

where the perturbation h^j,^ is a function of space and time. It will be easier for us to 
deal with the (spatial) Fourier components of h^^, and to avoid cluttered notation, we will 
subsequently denote the Fourier components h^i,{}i) simply by h^^. We choose to work in 
synchronous gauge, so /iqo = ^oi = ^io = 0, and the hij can be expanded in tensor harmonics, 
which satisfy V^Q^'") = -fc^QM |2|j3|], 

= Y: 'iHt^Q^'^H, + 24"^)glf , (13) 

m 

50 = Ex^™^Q^'"\ (14) 

m 

where Q*^™) and Qj-]"^ are scalar and tensor harmonics, respectively, and m denotes the 
"angular momentum" of the perturbation. For simplicity, we will write 

H^L^ = hi, H^T^ = Ht, 

H^T^ = hv, Hl = H. (15) 

For models with only scalar and tensor modes, /ly = 0. Our and are simply related 
to the variables used in Ref. |30] by 

h = QhL, r] = ~{hL + hT/3). (16) 

The perturbed stress energy tensor can also be broken up into scalar, vector, and tensor 
parts. Let us denote a cosmic fluid component (e.g., baryons, neutrinos, photons, cold dark 
matter, etc.) by index /. We then know that the stress-energy perturbations ST^-' are related 
to the perturbations 6p and 6p in the density and pressure, respectively, and to the velocities 
Vf and anisotropic stress ttj (see, e.g., Refs. ^'^^ more details), 

sT',=Y.j:iPf+Pf>f'Qi"'\ 

ST) = E E Sp^r^l'.Q^""^ + Pf^T^Qf^'. (17) 

If we consider only scalar and tensor perturbations, then the perturbed Einstein and Brans- 
Dicke equations are 

x' + 2-x' + k\ + 3/^10' = 1^ E i^Pf - ^spf) ; (18) 

(3 Z,UJ + O J 
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a (j) 



+ 



(19) 



/ 



(20) 



, „ a' , , 9 1 , X bvra 
/i^ + 2-/i^ - k^{hL + -/ir) = — - 
a 3 



pT^f — hrp \- k —; 



(21) 



a' . o SnGa'^ (2) 



a 



(22) 



D. Temperature and Polarization Anisotropics 



With these equations, one can find the evolution of perturbations using standard cos- 



mological perturbation theory; see e.g. Ref. p9|. The calculation of the CMB anisotropy 



runs in parallel to the one in the standard model detailed, e.g., in Refs. [pO| -p^. Here we 
summarize the procedure for such calculations. 

A temperature map T(n) of the sky (as a function of position n on the sky) can be 
expanded in spherical harmonics. 



Tin) 
To 



n 



(23) 



Im 



where the mode amplitudes are given by the inverse spherical-harmonic transform. Similarly, 
if we measure the Stokes parameters Q{n) and U{h) as a function of position on the sky, 
they can be assembled into a symmetric trace- free (STF) 2x2 tensor ||34|| , 



Ij Q(n) -U{n)sm9 
2 V -U{h) sine -Q{n)sm'^e 



(24) 



which can then be be expanded |Q, 
Vabin) 



(Im) (lm)ab 



(25) 



where the tensor spherical harmonics Y^^^^f^ and Y^^^^f^ form a complete basis for the gradient 



{lm)ab 



and curl components of the tensor field, respectively, and the multipole coefficients, 
and ct^J^) can be obtained by inverse transforms. 



(Im) 
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Thus, a combined temperature/polarization map is specified completely by the three sets 
of coefficients, aji^^, a^^^^, and a^j^y The two-point statistics of the T/P map are specified 
completely by the six power spectra Cf^'^' defined by 

(^flm)) = Sll'Smm', (26) 

for X, X' = {T, G, C}, and the angle brackets denote an ensemble average. Parity invari- 
ance demands that C'^'" = Cp*" = 0. Therefore the statistics of the CMB temperature- 
polarization map are completely specified by the four sets of moments, C^'^, Cj'~^ , C*P*^, and 
Cp*". These correlation functions are related to the ones used by Seljak and Zaldarriaga ||3^ 
by Cp^ = Cei/2, Cpc = Cci/2, Cj^ = Cci/V^, and our Cf^ is the same as their Cn- 

We can calculate the C;'s by convolving the initial metric perturbation power spectrum 

with the photon transfer function Ai(fc, Tq), 

Cj = {Any J edkP4k)[ATi{k)]\ (27) 
Cp^ = {An)' J edkP4k)[AGi{k)]', (28) 
= {AtcYJ edkP4k)[ATiA'i^!]. (29) 

The photon transfer functions Axi{k, tq) are obtained by integrating along the line of sight 
0, 



ATi{k, To) = / drS^k, T)ji[k{To - r)], (30) 
Jo 

^Gl{k,T,) = ^ £"dTSG{k,T)ji[k{ro-T)], (31) 

where ji{x) is the spherical Bessel function and ST,G{k, t) are the source functions describing 
the Thomson scattering of photons along the path, and 



5^ = ,(A^o + 2«' + | + H + ^)+e-V + a") 



^ ,(v, 3U'\ 3g"U 
_ 3g{T)U{T, k) 

- %{kTY ' ^^^^ 
n = At2 + Ag2 + Ago, 

where x = k{To — r) and a = {h' + Qri')/2k'. The visibility function g{T) is given by 
g = e""/?', where k{t) is the optical depth from conformal time r to the current epoch. In 
the Brans-Dicke theory, the derivatives of a are given by 
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fa'\ 








[a) 



a" = -2l-) a-2{-]a' + ri'^±^ a^p + p)a 




For the initial conditions on the scalar-field perturbation, we consider only the simplest 
case with ximt = Xinit — 0- The initial conditions for the matter are the same as those 



in the GR case As perturbations in the metric grow, perturbations in the Brans- 

Dicke field will also grow as shown in Eq. (|18|). However, for the initial condition we have 
chosen, the Brans-Dicke perturbation is so small that it has little effect on CMB anisotropy. 
An alternative choice of initial conditions for the Brans-Dicke perturbation would probably 
yield the same numerical results, because any initial perturbations are damped during the 
radiation dominated era [c.f. Eq. (|18D ]. 

The numerical calculation is essentially carried out by replacing the general-relativistic 
perturbation equations in a publicly available code [^,0 by those in Eqs. (p!9D -(p2D and 
including the evolution equation (^) for the Brans-Dicke field. In practice, there are a 
number of numerical issues and subtleties that arise, and some of these are detailed in the 
Appendix. 

We have chosen to work in the Jordan frame in which the equations for spacetime- 
metric perturbations are altered while the equations for the stress-energy perturbations are 
unchanged. We considered working in the Einstein frame, in which the metric-perturbation 
equations are unchanged, but found that the changes in the equations for the stress-energy 
tensor would be more difficult to implement numerically. 



III. BACKGROUND COSMOLOGY 

Let us now consider the background cosmology, and the boundary conditions for the 
homogeneous component of the Brans-Dicke field $ and its conformal-time derivative 
We define the cosmic scale factor at the present epoch to be = 1. In general-relativistic 
cosmology, the initial condition for the scale factor is a(r = 0) = 0. The conformal age of 
the Universe can be obtained by integrating, 

dr = da/a' . (36) 



JO 

For Brans-Dicke cosmology, additional boundary conditions are required for and 0'. 
One of these is determined by the requirement that the gravitational constant be in agree- 
ment with that measured today. This fixes 0] 

(37) 

at the present epoch. 

The cosmological solutions for Brans-Dicke theory have been studied extensively 
p0| , |37| -p9| . The Brans-Dicke field has a stiff equation of state; it dominates the dynam- 
ics at early stages of the expansion. However, for the era which affects the CMB anisotropy, 
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the Brans-Dicke field must be sub dominant, or else the expansion rate at nucleosynthesis 
would have been very different. Therefore, for a qualitative understanding of the expansion, 
we can assume that the change in does not affect the dynamics, and estimate how varies 
by assuming the Universe expands as in the GR case. 




X 

FIG. 1. The function Ki{x)/x versus x = m/T (increasing time). 

The equation of motion for the field is given by Eq. (^). It is analogous to a damped 
oscillator with a variable friction force. The initial "velocity" 0' is damped in a few Hubble 
times. Therefore, for most of the time concerned, would only vary slowly. The right-hand 
side of Eq. ([T0| ) is the "driving force" for the motion of 0. It is proportional to p—3p. If p—3p 
were to vanish during radiation domination, then the quantity y = a?(f)' would be constant. 
However, it does not vanish during the radiation dominated era; in fact it is greater, even 
though the ratio (p — 3p) /p is smaller. There are two kinds of contributions to p — 3p. First, 
the non-relativistic matter, including both baryons and cold dark matter, always contributes. 
Second, if one massive relativistic particle is present, it always dominates.^ For one species 
of massive relativistic particles. 



^Here we do not consider the very early era during which everything, including baryons and 
cold-dark- matter particles were still relativistic. 
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where f{p) is the distribution function. For example, for a Boltzmann gas with zero chemical 
potential, f{p) = e~^^^, and the result can be expressed in modified Bessel functions, 

g Kiim/T) , , 

p-3p= , , L ■ (39) 

Fig. plots this function. As one can see, as T — > oo, m/T 0, and this function rises 
rapidly. 

In the present paper, we will not consider massive neutrinos. In a CDM model, the last 
decoupled massive relativistic particles are electrons and positrons. They annihilate below 
T ^ rrie = 0.511 MeV. For the scale we are interested in, the main contribution comes from 
the cold dark matter, which scales as pc = pcoa~^. 

After e~^e~ annihilation, in the radiation dominated era, 

a Qc T, (j)' ^ ci + C2a~^, (40) 

and one can see that approaches a "terminal velocity" and the initial velocity quickly dies 
out. From T ~ 0.5 MeV to T ~ 10 eV (matter-radiation equality), the initial velocity of </> 
is suppressed by a factor of 10^, and this initial "velocity" is constrained by nucleosynthesis, 
so it cannot be too large. We estimate, 

2u + 3 fcb'V ( a'^ ^ 



at the end of nucleosynthesis. We find that for all practical purposes, we can take 

= 0. (42) 

This is effectively the Brans-Dicke initial condition 0] proposed in their first paper. 

In the matter-dominated era, (p varies as oc a^^^^'^^\ For models with a; > 0, the value 
of (f) increases with time, whereas for models with u < 0, (p decreases with time. Fig. |^ 
shows the evolution of (p. 

We also note that in Brans-Dicke theory, the matter density is not precisely equal to the 
critical density in a fiat Universe. The critical density in the Brans-Dicke theory depends 
on the parameter u. If we still define our relative densities in the usual way, i.e., 

= Pi/pc, Pc = 3H^/S7iG, (43) 

then we will have fl = ^ fo^ Acit geometry. Let us define 

; (44) 



4> 



a=l 



then from Eq. (IS) we have 
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With the matching condition, cf) = {2uj + A)/{2uj + 3), one can obtain the value of f2 if D is 
known. To proceed, one may start with some value of p and solve the evolution equation to 
obtain Hq and D. In practice, for the models considered here the difference is very small. 
We have tested for a few models with u equal to a few hundred, and it does not make any 
significant difference. 

IV. RESULTS 

We now illustrate our numerical results with a few representative models. First we 
consider a COBE-normalized flat standard CDM (SCDM) model, with n^, = 0.03, = 0.97 
and h = 0.65. The angular power spectra for u = 200, 500, and —200 are plotted in Fig. 
^ For comparison, we have also plotted the general-relativity result with the same physical 
parameters in the same plot. 

For uj = ±200 the difference between Brans-Dicke models and general relativity are 
clearly discernible As can be seen, both the normalization and width of the acoustic peaks 
are changed. The Brans-Dicke model with a positive u has higher and broader acoustic 
peaks, while the negative-ci; model has lower and narrower peaks. We have checked that the 
perturbations in the Brans-Dicke field near the time of decoupling in this model are very 
small. Thus, change in the acoustic-peak structure is due primarily to the change in the 
expansion rate near decoupling. For u = 500 the difference is much less pronounced. The 
polarization spectra are similarly affected. 

The Brans-Dicke field also affects the transfer function. Fig. ^ compares the matter 
transfer function in a general-relativity model, a Brans-Dicke model with u = 200, and one 
with UJ = 50. For the Brans-Dicke models, the bend of the matter power spectrum occurs at 
shorter wavelengths, and there is thus more small-scale power, in agreement with the claims 
of Ref. 

The CMB power spectra are also affected by other cosmological parameters, and it is 
possible that variation of some other parameters might mimic the effect Brans-Dicke gravity. 
For example, we plot the Brans-Dicke model with u = 200 and Qb = 0.030 along with a 
general-relativity model with Uf, = 0.032 in Fig. ^. This general-relativity model mimics 
the Brans-Dicke model up to the first acoustic peak. Note, however, that this different flf, 
model does not fit the polarization better — in fact, the fit for the polarization is even worse. 
Therefore, observation of the polarization may help to lift this degeneracy of parameters. 

It may also be possible to mimic the Brans-Dicke model entirely with a general-relativity 
model by adjusting more than one parameter. To investigate properly the possible degener- 
acy of the effect of varying uj with the possible effect of varying some combination of other 
cosmological parameters, we calculate the covariance matrix This also allows us to 

simultaneously estimate the precision with which u (actually uj~^) can be recovered with 
a CMB map. We first consider only a temperature map and later consider the additional 
information that comes from the polarization. If the true parameters which describe the 
Universe are given by Sq, then the Fisher information matrix is defined by 
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FIG. 3. CMB temperature/polarization power spectra for flat SCDM models in general rela- 
tivity and in Brans-Dicke theories with lo = 200, 500, and —200. 
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FIG. 4. Spatial power spectra for SCDM models in Brans-Dicke theory with 
= 50 and for general relativity. 
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FIG. 5. CMB power spectra for a Brans-Dicke SCDM model with fi;, = 0.03 and oj = 200, and 
for general-relativistic models with Vti, = 0.3 and = 0.032. 
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If the observed Q's are nearly Gaussian distributed around Q(so) with variance ai, the 
covariance matrix [C] = [a]~^ gives an estimate of the standard errors that would be ob- 
tained from a maximum-likelihood fit to the data. Approximately, the standard error with 
which the parameter Sj could be obtained (after marginalizing over all other undetermined 
parameters) would be cxs- = C}/'^ . 

Consider a CMB experiment that maps the temperature of the entire sky with a Gaussian 
beam of width ^fwhm, with a noise per pixel of (jpix. If a fraction /sty of the sky is used after 
a foreground cut, and the noise in each pixel is uncorrelated, then the the standard error 
with which each C',^^ can be recovered is 

1/2 

(G + ^^V'^^'), (47) 

where cxf, = 7.42 x 10~^(6'fwhm/l°), and the inverse weight w^^ is given by 

= 4vraJJiVpi„ (48) 

where cXpix is the noise per pixel, and iVpix — 40, 000(6'fwhm)~^ is the number of pixels. 

The goal of the MAP mission is to measure the temperature anisotropy with 6'fwhm = 0.3° 
and (Tpix = 20 /xK, which corresponds to (^assummg a one-year experi- 

ment). The Planck Surveyor has a mission goal of 6'f„hm = 0.1° and CTpix = 5/iK, which 
corresponds to = 6 x 10"^^. Assuming /sky = 0.67, we calculate the Fisher infor- 
mation matrix with ^fwhm = 0.3° and 0.1° for a variety of values. The results are 
shown in Figs. ^ and ^. Our fiducial model is a COBE-normalized flat CDM model with 
{ho, Qb, ns, T, l/uj, {Q / QcobeY} = {0.65, 0.03, 1, 0.5, 0, 1} where ns is the primordial power- 
spectrum index, r the reionization optical depth, and Qcobe('^s = 1) = 18 /iK is the COBE 
normalization We assume 3 generations of massless neutrinos and no massive neutri- 
nos, and consider scalar modes only. The derivatives of Ci are calculated by varying each 
of the parameters by 0.5%. To calculate the partial derivatives of the Q's with respect to 
1/a;, we compare the general-relativity model with a Brans-Dicke model with u = 200. Our 
calculation sums up modes up to / < 3000. We have checked to make sure that our results 
are not sensitive to the step size for calculating the derivative nor the cutoff of I. 

The value of cri/t^ depends on the number and uncertainty of other parameters. Table 
I lists the standard errors that could be obtained by marginalizing over all others for the 
various parameters we consider. Fig. || plots the smallest value of u that could be distin- 
guished from u = CO (i.e., general relativity) at the 95% confidence level. So, for example, 
if all the parameters listed above were unknown and had to be determined from CMB data 
alone, then the CMB would be marginally competitive with (current) solar-system tests; 
i.e., if we marginalize over Qcobe, ^b, h, n^, and r, then the smallest value of u that could 
be distinguished from oo is ~ 100 for MAP and ~ 800 for the Planck Surveyor. 

On the other hand, if we assume that all parameters except for u and the normalization 
can be determined completely from other experiments, then the sensitivity to a finite uj can 
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FIG. 6. The largest finite value of iv that could be distinguished from infinity (i.e., general 
relativity) at the 2a level as a function of the pixel noise of a given experiment that covers two-thirds 
of the sky. The fiducial model here is a standard CDM model (no cosmological constant). We show 
results for two beamwidths, ^fwhm = 0.1° and ^fwhm = 0.3°. We assume here that Us, r, and 

Q are marginalized over. The solid curve corresponds to ^fwhm = 0.3° with temperature data only, 
and the short dashed curve also includes the polarization. The long dashed curve corresponds to 
^fwhm = 0.1°with temperature data only, and the dotted curve includes also the polarization. The 
expected values of cjpix for MAP and Planck are indicated. 





ho 




ns 


T 


(Q/Qcobe)^ 


1/uJ 


value 


0.65 


0.03 


1.0 


0.5 


1 





'7MAP(7') 


0.045 


0.0054 


0.043 


0.057 


0.057 


0.0050 


CTMApiT + P) 


0.031 


0.0036 


0.031 


0.035 


0.053 


0.0034 


<y Planck 


0.0045 


0.00049 


0.0081 


0.013 


0.018 


0.00062 


0'Planck(r + P) 


0.0037 


0.00040 


0.0055 


0.006 


0.017 


0.00049 



TABLE I. Error estimates for parameters of an SCDM model. Here, MAP is assumed to have 
^fwhm = 0.3° and w^^ = 2 x 10^^^; Planck is assumed to have ^fwhm = 0.1° and = 6.3 x 10^^'^. 
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FIG. 7. Same as Fig. ^, except that all other parameters (except the normalization Q) are 
assumed to be known. 
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be improved, as also illustrated in Fig. |^. For example, a finite value of uj as large as ~ 500 
could be detectable with MAP and ~ 2500 for Planck. 

The future satellite missions will measure not only the temperature anisotropy but also 
the polarization. It is possible to improve the accuracy of cosmological-parameter deter- 
mination by combining the temperature and polarization data, as Fig. ^ illustrated that 
polarization may help break degeneracies in parameter space. 

To include the polarization data, we generalize Eq. (EBI) to 



dC 



dsi 



^-1 



da 



Y,l 



(49) 



X,Y I I ^'^^ '^'^J J 

Here, X, Y =TT, GG, CC, TG, and [S^^]x.y are elements of the inverse noise covariance 



matrix S. The elements of S were given in Refs. [0,^]. If the two linear-polarization states 
are given equal integration times, the total number of photons available for temperature 
measurement are twice of that for polarization measurement, thus 



= 7^K.f- (50) 



If the number of pixels are equal, then 

w^^ = Wp^. (51) 

The results obtained from combining the temperature and polarization data are also 
plotted in Figs. |^ and |^, and these plots show that the sensitivity to the Brans-Dicke pa- 
rameter u could be further improved by including polarization. The effect is particularly 
strong when the effect of Brans-Dicke gravity is degenerate with that of variation of other 
cosmological parameters in a temperature map. By including the polarization data in MAP, 
the CMB should be sensitive to models with uj < 150 at 95% CL when all other parameters 
are undetermined, ot uj < 500 when only CMB normalization is undetermined (in this case, 
there is not much gain from polarization). For Planck, these numbers are u < 1000 and 
UJ < 3200, respectively. On the other hand, if all the other parameters are known, the benefit 
gained from adding polarization measurement is less obvious; only for detectors with pixel 
noise less than 5/iK there is a difference. 

There has been much recent interest in models with a nonzero cosmological constant 
prompted, in particular, by the evidence for an accelerating universe from supernovae 



We have also performed our calculations for a ACDM model, and the results are shown in 
Fig. H for u = 200. The results for the Fisher-matrix analysis are shown in Figs. ^ and |l^, 
and they are similar to those in the SCDM case. 

Finally, we have noted before, if uj is finite, then Q is not precisely equal to unity [cf., 
Eq. (P^l)], and one might wonder whether the effects on the CMB power spectra of varying 
u can be mimicked by varying Q. We have checked that for the Brans-Dicke models we have 
investigated, the change in the CMB power spectra is much too large to be attributed to 
this shift in Q. 
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1 

FIG. 8. CMB power spectra for ACDM models in Brans-Dicke theory with to = 200 and in 
general relativity. 
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FIG. 9. The same as Fig. |6|, except that here we use a flat cosmological-constant model with 
a nonrelativistic-matter density 0,q = 0.4. 
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FIG. 10. Same as Fig. except that all other parameters (except the normalization Q) are 
assumed to be known. 
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V. CONCLUSION AND SUMMARY 



We have developed a formalism for calculating the CMB anisotropy in cosmological 
models with Brans-Dicke gravity. This was done by modifying standard Boltzmann codes 
for CMB power spectra. Because Brans-Dicke theory satisfies the medium-strong equivalence 
principle, only those equations determining the evolution of the metric need to be modified; 
the equations of motion for various matter or radiation components are the same as in 
general relativity. One boundary condition for the Brans-Dicke theory is determined by 
requiring (p = at the current epoch. Another is given by adopting the Brans-Dicke 
initial condition a^0' = at early time (after the annihilation of electron-positron pairs). 

This formalism is then used to calculate the CMB power spectra in several models. 
We find that in Brans-Dicke models, both the height and width of the acoustic peaks are 
changed. While there is some degeneracy with different cosmological parameters at low / in 
a temperature map, we demonstrate that the effect can be distinguished by going to higher 
acoustic peaks and by observing the polarization of the CMB. Our results show that with 
high-quality CMB data, the CMB anisotropy may provide a powerful test for Brans-Dicke 
theory that is competitive (and complementary) to solar-system tests. 

As an example, we examined a flat SCDM model. MAP temperature data should be 
able to distinguish Brans-Dicke gravity with u; < 100 from general relativity at the 95% CL 
if all other parameters must be simultaneously determined from the CMB, or a; < 500 if all 
other parameters except for the CMB normalization are fixed. With Planck, these numbers 
are 800 and 2500 respectively. 

Furthermore, even better results are achievable if both temperature and polarization 
data are used. For MAP, the two limits are raised to 150 and 800, respectively, and for 
Planck, 1000 and 3200, respectively. We also examined the case of a fiat ACDM model and 
found similar results. 

In conclusion, the differences between the CMB power spectra expected in general rel- 
ativity and those in Brans-Dicke models with acceptable values of uj are small. However, 
our Fisher-matrix analysis shows that if systematic effects can be controlled, then the CMB 
sensitivity (from Planck) to a finite value of cu might be competitive with that from solar- 
system tests. We re-emphasize that the CMB will provide a new and independent test of 
gravity in stronger fields and at earlier times. Thus, it is conceivable that the CMB will 
provide a unique test of some scalar-tensor theories in which uu would have been smaller at 
earlier times. 
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APPENDIX A: NUMERICAL IMPLEMENTATION 



In this Appendix we briefly describe the numerical implementation of the calculation. 
First we consider the background evolution. The boundary conditions for the Brans-Dicke 
field (f) are given in Eqs. (|37|) and (^21) . From Eq. (|37D , we know the end-point value of 
but not the initial value. We pick an arbitrary epoch with temperature 10 eV<^ T ^ 0.5 
MeV, then pick some value of with 0' = and integrate forward until the scale factor 
a = 1.0 This process is reiterated with different trial initial values of until the condition 
00 = (2^1^ + 4)/(2a; + 3) is satisfied to the required precision. Formally, the process is 
equivalent to numerical root finding, and we use a Brent algorithm |Q to find the root. 



In the original CMBFAST code, any epoch in the evolution of the Universe is specified by 
the cosmic scale factor a, and the time r is obtained by integrating Eq. (PB]). In Brans-Dicke 
theory, the value of a' /a is given by Eq. (^, and it is no longer convenient to use a as the 
argument. Instead, we use r to specify the cosmic time, and the scale factor a is obtained 
by solving the whole set of background-evolution differential equations. This is done in the 
beginning of the calculation, and the values of {T,a{T)^a'{T),(l){T),(f)'{T)} are then stored. 
Subsequently, given either r or a, the whole set of these values corresponding to that epoch 
can be obtained by lookup and/or interpolation. The second-order derivatives a" and 0" 
can also be obtained by numerical difference. 

This more complicated implementation of the cosmic history also demands modification 
in other parts of the code. In the original code, the expansion rate is calculated from the 
density of the universe using the Friedman equation; in this new code, we replace it by 
interpolation from the stored data in all such cases. 

For example, in the original code, the baryon temperature, ionization fraction, and 
baryon sound speed are given as a functions of a, and are calculated by a simple inte- 
gration of the Friedmann equation. In the new scheme, this calculation is modified so that 
the value of a' is obtained by lookup or/and interpolation of the stored data. We also modi- 
fied the code for the calculation of the baryon sound speed. In the original calculation, there 
are occasional jumps in the baryon sound speed, which may be due to the truncation error 
in calculating the derivative dTf,/dlna. We have modified the algorithm so that there is no 
jump in this calculation. However, our tests show that these occasional jumps do not have 
any significant effect on the end result, probably because only interpolated values are used, 
and the result is mostly important only in a limited range. 

For the calculation of the perturbative part, we note that in Brans-Dicke theory, the 
equation of motion for the matter or radiation are the same as in general relativity. All we 
need to do is to replace the perturbed Einstein equations by Eqs. (|I8|)-(^2D. Following the 
original code, we use Eq. (|T^ to force conservation of energy and reduce the numerical error 
in solving the ordinary differential equations. We have tested that when ^ oo we recover 
the general-relativistic result produced by the standard code. 

Finally, the temperature and polarization anisotropy may be obtained by integrating 



^In principle, one may also integrate backwards from the present epoch to this early stage and 



thus obtain the "initial value", but this procedure is susceptible to numerical instability |17]. 
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Eqs. (p7|)-p9|). The expression for the source function is the same, but note that there are 
Brans-Dicke corrections to the derivatives of metric perturbations as given in Eqs. (|3^ - (^5|) , 
and these must be implemented in the code. 
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